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Abstract. We study the coupled rotational diffusion in a two-particle chain on the basis 
of a Smoluchowski equation and calculate time-correlation functions that are measurable in 
an experiment. This might be used to explore hydrodynamic interactions in the limit where 
lubrication theory is valid. 
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1. Introduction 

Colloids are widely considered as models for atomic systems (l] |2l . However, a specific 
signature of colloidal suspensions are hydrodynamic interactions (3j |4l . Moving colloids 
interact with each other through the flow fields they create. This is a true multi-body problem 
which can only be handled by approximate methods such as multipole expansions for large 
particle distances (see, e.g., Refs. |4l|5|) and lubrication theory when they come close to each 
other f6l. Conventionally, hydrodynamic interactions are monitored through their effect on 
self and collective diffusion in colloidal suspensions 01 IH ^ but recent experiments with 
optical tweezers on a pair of particles |8l also allow a controlled exploration of hydrodynamic 
interactions as a function of particle separation confirming standard approaches due to Oseen 
and Rotne-Prager 1 3 , 4 1 . 

Recent work also studied the rotational diffusion of tracer particles 191 1101 [TTl or colloids 
trapped in optical tweezers 1 12|. An experimental system introduced by Bibette et al. 1 13 1 
suggests a possibility to directly measure the effect of hydrodynamic interactions on the 
rotational diffusion in the limit where lubrication theory is valid. Charged superparamagnetic 
particles under the influence of a magnetic field form chains where the particle separation 
can precisely be tuned by the magnetic field strength. By labeling the particles with a 
phosphorescent dye flTll or by using birefringent colloids fl41 . rotational diffusion can then 
be monitored. 

The present article investigates rotational diffusion in a two-particle chain theoretically. 
It first introduces the Smoluchowski equation to treat rotational diffusion and then discusses 
observabels to be measured in an experiment. 



2. Rotational Diffusion and the Smoluchowski Equation 



Let us first shortly review the rotational diffusion of one particle that can also be considered 
as a random walk on the unit sphere. The probabiUty density P{C', t) of finding the particle 
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with an orientation given by the unit vector i> at time t satisfies the Smoluchowski equation 

m 

'-^=D,VlPiO,t) . (1) 

This is in complete analogy with translational diffusion. However, instead of the Laplace 
operator, the square of the nabla operator in angular space has to be used, where the index 
r means rotation. Since Vr = where L is the angular momentum operator known 

from quantum mechanics, all the algebra developed for L |15| is also valid for V,-. To be 
concrete, we note that 



d 



^^dO^\ I . (2) 

dip 

where 0, -d are the spherical coordinates to represent z>. The rotational diffusion constant Do 
in Eq. Q is related via an Einstein relation to the mobility /io = l/{8TTr]a?) for rotational 
motion, i.e., Z?o = kBTiiQ. Now, since the operator possesses the spherical harmonics 
y,'j(i>) as its eigenfunctions, V^y^j(/>) = + l)y^(i>), one calculates, with the same 
method as presented below, the time-correlation functions 

{Y::[0{t)]Yl,MO)]) ^ . (3) 

Here C'{t) means orientation of the particle at time t and the symbol * means complex 
conjugate. For small times, one shows with the help of Eq. Q (for details, see below) that 
the square of the angular displacement of /> exhibits the typical diffusive behavior: 

{\u{t)-i>{0)\^) ^4Dot . (4) 

We now concentrate on a chain of two particles of radius a whose centers are connected 
by the vector r. For simplicity, we only consider their rotational diffusion, i.e., we disregard 
any coupling to translational motion. We especially take the particles' separation r as fixed, 
e.g., by assuming that fluctuations around the equilibrium separation, governed by the two- 
particle potential, are negligibly small. The coupled rotational diffusion is described by self- 
diffusion tensors Dn ~ D22 for particle 1 and 2, respectively, and the tensors D12 — D21 
encoding the hydrodynamic interactions between particle 1 and 2. As in the one-particle 
case, these quantities are related to mobilities by an Einstein relation, Dij = ksT^^y The 
mobilites ^.^ connect the torque Tj on particle j with the angular velocity uii of particle 
i: LJi — fiijTj. Due to Lorentz's reciprocal theorem |3|, they fulfill fl^j ~ where t 
means transposed matrix. If, in addition, the particles are identical, fi^^ — fi^^. The uniaxial 
symmetry of the two-particle chain determines the form of Dij: 

= A^l + AD.yf (g) r with AL>,y ^ ~ D^^ , (5) 

where 1 means unit tensor, f = r/r, and ® means tensor product. The constants d\- and 
-D^ refer, respectively, to rotational diffusion about the two-particle axis f or a direction 
perpendicular to it. The related mobilities as a function of reduced particle distance r/a are 
plotted in Fig. [2 The library HydroLib 1 16 1 allows to calculate their values ranging from 
small particle distances where lubrication theory has to be applied to large distances where 
expansions into a/r are applicable. For comparison, the Rotne-Prager approximation is also 
shown. Both graphs of Fig. ^ illustrate that the Rotne-Prager approximation works well for 
center-to-center distances larger than 3a. To see noticable deviations of the mobilities /i^^" 
and jiYi' from the single-particle value — l/{ST:-qa^) (the index rr refers directly to 
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Figure 1. The rotational mobilities of the two-particle chain in units of fig = 1/ (SvrTja'^ ) as a 
function of reduced center-to-center distance r/a. Full lines: values calculated with numerical 
library HydroLib; dotted lines: the Rotne-Prager approximation. 



the rotational degree of freedom), the particles have to be close. Note that the mobilities in 
both graphs stay finite when the particles approach contact at r = 2a; the derivatives of the 
perpendicular coefficients, however, are singular. 

The Smoluchowski equation determines the temporal evolution of the probability density 
P(i>i, i>2, t) to find particles 1 and 2 with respective orientations given by i>i and i>2 at time 

ot 

where 

is = Vrl ■ DnVrl + Vr2 ' £>22V^2 + 2Vrt • -Dl2V^2 (7) 

denotes the Smoluchowski operator that is a generalization of -Do in the single-particle 
equation 0. Since both particles are identical and with the help of Eq. (|5j, it is rewritten as 



d d 

+ 2Di^Vrl-Vr2 + 2ADi2 ^ ^ (8) 

0(pi 0(p2 

where d/dtpi is the z component of V,i with z\\r. In analogy to the wave function determined 
from Schrodinger's equation in quantum mechanics, the time evolution of the probability 
density P{C>i,C>2,t) is known in principle when eigenvalues and eigenfunctions of the 
Smoluchowski operator are known. Due to the analogy with the angular-momentum algebra, 
the eigenvectors of the unperturbed problem, i.e., particles not coupled by hydrodynamic 
interactions {D^2 = AD12 — 0), are just products of two spherical harmonics: <l>(i>i, 1^2) = 
('^i)^m2 ('^s)- They are even eigenfunctions of the last term in the second line of Eq. (|8}. 
The operator V^i • ^r2, however, mixes the eigenfunctions of the unperturbed problem. How 
this is done, can be calculated with the help of "ladder operators" introduced in full analogy 
to the angular momentum algebra 1 15J: 

= V„ ± iVry (9) 

AppHed to spherical harmonics, they yield 



VfYiiu) = i^l{l + l)-m{m±l)Yi^,{0) = *qtll±i(i>) . (10) 
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We rewrite the crucial term V^i • Vr2 of Ls with the help of 

Vrx = ^(V++V-) and Vrj, = ^(V+- V7) (11) 

as 

V.i . v.. = i(V+ V;, + V;,V+ ) + (12) 

and can now fully determine how the Smoluchowski operator acts on the unperturbed 
eigenfunction: 

- LsY^J^^ = {Di^Mh + 1) + Hh + 1)] + Ai^nK + ml) + 2D\^m^m2}Y^J^^ 
where we used = + AI?i2 and the coefficients are defined in Eq. ( fTot . 



3. Time Correlation Functions and Observables 



We are interested in quantities that can be measured in experiments. We therefore define the 
most general time correlation function 

= (>;[^:[^>iW]i;j?;[^2W] y^, [z>i(o)]y;4[«>2(o)]) (i4) 

which can be related to useful observables. Formally, F is calculated using the propagator 
P(i>i, £>2, t\i>i, O2, 0) that gives the probability of finding the two particles with orientations 
i/i and C>2 at time t when they had with certainty the orientations 0[ and O2 at time t = 0: 

r(;S:?|t)= / jY^:{O,)YX*{O2)P{i>,,i>2,t\O[,O'2,0) (15) 

X T4^(i>'i, i>^, t - 0) Y^, (i>'i) rj, (i>^) dO,d02dO[du'2 

where W{i''i,i'2,t = 0) is the probability distribution for £>[ and i>2 ^t t = 0. In the 
following, we will use an isotropic distribution W = l/(47r)^ since the interaction potential 
of the particles does not depend on their orientations. However, one could also think about a 
situation where one first "aligns" the particles and then lets them evolve with time t. In such a 
case, W would be given by a product of delta functions. The time evolution of the correlation 
function is calculated from a master equation that we derive by taking the time derivative 
of Eq. il5i . then using the Smoluchowski equation ^ for the propagator and finally letting 
Ls — Lg act on the spherical harmonics at time t: 

|r(r:;?:'::?i^) = / j[LsY^:{ui)Y^:{u2)]p{o,.02Ai^\,u'2,^) 

xW{0[,0'2,t^O)Y'},y,)Y^,^{0'2)dO,du2du[du'2 . (16) 
With the help of Eq. ( I13l l and the defintion ( I15t . the master equation assumes the form 

- |r(^};?™;:?|t) ^ {Dt.Uh + 1) + hih + 1)] + ADniml + ml) 

+ 2Dlm^m2}TC^^rZl\t) (17) 

4- n-L [^+ r-^ -n(hl2mi+lm2-l\,\, - + T^(hl2mi-lm2 + l \,\] 
"I ^12Rimi^/2m2-^ W'll^mim^ ^ '-hmi i2m2 Wi/^mim^ r^J 

The time evolution of various observables can now be calculated with the help of this equation. 
We illustrate two cases which should be measurable in experiments. 
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Figure 2. The effective rotational diffusion constants D'^, D^^ + dIj/S, and - 
in units of Do as a function of the reduced particle separation r/a. 



3.1. One-particle diffusion 

Let us investigate the time correlation function for the orientation of particle 1 : 

m— — 1 

The first equality is just the addition theorem for spherical harmonics 1 15| and the second 
equality used h = ^2 = — ^^'2 — 0' i-^-' ^0" — l/V^ for the second particle in the 
definition (O of r(. . . \t). Since c+''" = 0, r(}[S;;^[;|t) does not couple to other T's. The 
evolution equation is therefore simple and gives 

r(S|i) = r(}Slo)e-(^^"+^^"'"^)* (19) 

Compared to the single-particle result of Eq. ^ {I ~ 1), the decay rate of the correlation 
function now also depends on the azimuthal quantum number m due to presence of the second 
particle. With r(}Q™Q|0) = l/(47r)^ [isotropic distribution of i>i(0)], the correlation function 
fWi becomes 

(i>i(t)-i>i(0)) = ie-2^i"i'(l + 2e-^^"*) . (20) 

Finally, the mean-square angular displacement of Oi{t) is calculated from (|i>i(/:) — 
'^i(O)P) = 2(1 - (i'lit) ■ «>i(0))) and for small times it reads 

( \u,it) - t>i(0)p )^4Df,t with Df, = + . (21) 

So Dq in the analogous single-particle equation (0} is replaced by the average of the self- 
diffusion constants D^i and that encode the presence of the second particle. The effective 
diffusion constant _Djf in units of Dq is plotted in Fig.|2]as a function of the reduced particle 
separation r/a. 
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3.2. Two-particle diffusion 

To access the diffusion constant I^fj, we consider the correlation function 

mi— — 1 m2— — 1 

The correlation functions r(}}^^^^^|t) do not couple to other F's and therefore relax with 

a rate 2{Dii + D^-^^ + 0^2) as determined from the first two lines of Eq. ilH . The 
quantity Till°^^^Jt) couples to r(Ji^-J„Jt) and r(}}-ii„Jt) but can be determined 
straightforwardly from the three coupled evolution equations. So, ultimately the correlation 



function (I22> appears as a sum of three exponentials. Since we are mainly interested in 
the short-time limit of Eq. i22i . we present a shortcut towards the result. With the Taylor 
expansion 

-'^(llT?X2m2 10 ~ ^^jj.^2 '^™l"i2 ~ 0,mit , (23) 

where the coefficients ami ^re directly determined from the master equation ( I17> using 
r(n™-™2 10) - r(}J-ii„jO) = 0, we immediately arrrive at 

{O,{t)-O2{t)O,{Q)-O-2{0))-^[l-A{D<if + Dl/m ■ (24) 

Note that Df^, familiar from the one-particle diffusion [see Eq. (I21H . is modified here by 
dI^/S. Figure|2lillusti-ates Dl^ + dI^/S as a function of particle separation. 
In the same manner, we can also determine the short-time limit of 

{O,{t)-O,{0)O2{t)-O,{0)) = {fy E E ^Ci'ZZm (25) 

?n 1 — — 1 m2 — — 1 



and obtain 



with 



[Oiit) ■ i>2(0)i>2(<) • i>i(0)) « ^ [1 - 4{Df, + 01^)1] (26) 



^cff ^ . (27) 

We again plot Z3jf + Djf in Fig.|2] A comparison of all three effective diffusion constants in 
Fig.|2]only reveals small differences that pose a challenge for measurements in an experiment. 
Nevertheless, for small particle distances, the effect of hydrodynamic interactions calculated 
on the basis of lubrication theory should be measurable. 

A natural extension of the theory presented here is a chain of more than two particles. 
However, the mobilites or diffusion tensors in such a system cannot any longer be described 
by the two-particle system; when a particle is situated between two other colloids, three- 
body effects become important, certainly for small separations. So experiments could be 
used to measure deviations from the effective diffusion constants presented in this article. 
Furthermore, the coupling to positional fluctuations of the particles has to be incorporated. 



2Dt2 + dI 
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